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Abstract 

Results of a next-to-leading order calculation of the inclusive cross section for a photon and a 
heavy quark (charm or bottom), pp/pp — > 7 + Q + X , are presented. Pointlike photon subprocesses 
through 0(aa 2 s ) and fragmentation subprocesses through G(o? s ) are included. The calculation is 
performed using a phase space slicing technique so that the effects of experimental cuts can be 
included. Results for the ratios of the charm and bottom cross sections are presented and the 
systematics of the various subprocesses for both the Tevatron and the LHC are compared and 
contrasted. 
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I. INTRODUCTION 



Large momentum transfer processes have a long history of providing information on the 
substructure of hadrons, the nature of their constituents, and how they interact [1!] . Photons 
provide an excellent probe for such purposes due to their pointlike electromagnetic coupling 
to the quark constituents of the colliding hadrons. The production of large momentum 
transfer photons has played dual roles of providing information of parton distribution func- 
tions (PDFs) |2| and testing the adequacy of the perturbative techniques used to calculate 
the hard scattering subprocesses [j], |j, [l]] . 

The single photon cross section, either inclusive or subject to photon isolation cuts, 
provides the basic observable for direct photon studies. The calculation of this cross section 
involves integrations over the phase space variables of the accompanying partons, thereby 
limiting the information which can be obtained about the underlying subprocesses. More 
information can be obtained if one can measure an associated jet as has been done recently 
by the D0collaboration [6]. Even here, however, one is summing over many subprocesses 
with various flavors of partons. Additional information can be obtained if the flavor of the 
produced jets is tagged. Identifying jets which contain a heavy quark provides exactly this 
opportunity. 

In this paper we investigate in detail one particular piece of the direct photon calculation, 
namely the associated production of direct photons and heavy quarks, where the heavy 
quarks are either charm or bottom. Some of the contributing subprocesses are dependent on 
the charge of the heavy quark while others are not. By considering both charm and bottom 
quarks one can examine the relative roles of the two contributions. In some kinematic regions 
the results are dependent on the heavy quark PDFs, opening the possibility of testing the 
current calculation of such PDFs. 

New measurements of this process by the D0 and CDF groups are in progress and should 
be available in the near future. A comparison with these measurements will be presented in 
a forthcoming paper. 

The production of heavy quarks at high-p^ has the potential to generate logarithms of 
the form ln(pT/fnQ) as a result of collinear configurations involving Q — > Qg and g —* QQ. 
These logarithms can be resummed via the DGLAP equations for appropriately defined 
PDFs and fragmentation functions (FFs). This is commonly referred to as a variable flavor 
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scheme with either four or five flavors. In such schemes the heavy quarks are treated as 
massless. The dominant remaining mass effect is due to the imposition of a threshold 
constraint such that the PDFs and FFs are taken to be zero when the hard scattering scales 
are smaller than the quark mass. The calculation presented here has been performed using 
the variable flavor scheme. 

There have been previous studies of this process Q,y,y,Q. in Ref. 0,3 results are 

shown for the production of a direct photon plus charm. Here we also provide results for 
direct photon plus bottom production, and a comparison between the charm and bottom 
case. As noted above, this comparison depends on the relative roles of terms which depend on 
the heavy quark charge and those which do not. We also extend the calculation by treating 
the photon fragmentation contribution to next-to-leading-order (NLO). The previously cited 
references treated this component only in leading order (LO). One further technical point 
is related to the treatment of final state collinear singularities in the case when a gluon is 
emitted collinearly to a final state heavy quark. In Ref. [?| this singularity is factored into a 
charm FF. The present calculation is for the case of a photon produced in association with 
a jet which has been tagged as containing a heavy quark, so the fragmentation function is 
replaced by an appropriate jet definition. 

In order to be able to work in the massless approximation the heavy quarks and photons 
produced need to carry a transverse momentum, px, which is few times larger than the 
mass of the heavy quark tuq, i.e. pt > 10 GeV. Since the lower bounds for the values of 
the transverse momenta for direct photons and heavy quarks measurable at both the D0 
and CDF collaborations at Fermilab are above pt > 10 GeV, a comparison with a massless 
calculation is appropriate. If however we are close to the threshold region for production 
of heavy quarks, i.e. rriQ ~ pr, their mass needs to be retained, and they are assumed to 
be only produced externally in pairs, as end products of the hard-scattering. This is called 
the Fixed Flavor Number scheme (FFN), as the number of flavors that compose the nucleon 
remains fixed and it does not depend on the center of mass energy of the hard scattering. 
Here the proton is assumed to be composed only of light flavors, and in lowest order there are 
only two subprocesses in which the direct photon and heavy quarks can be produced. These 
are gg — > •yQQ, and qq — > jQQ. A study of the case when the quark masses are retained has 
been done at LO [9|, where a comparison between the LO massive and massless approaches 
has been shown to give very similar results. In Ref. B] a differential cross section in the 
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transverse momentum of the photon up to values of pr-y ~ 50 GeV is presented. There the 
difference between the two approaches in the LO is minimal. 

This paper proceeds as follows: in Section II a description of the theory and techniques 
for the calculation are outlined. In Section III results for the differential cross sections are 
shown. Predictions for both the Tevatron and LHC are presented and compared. The 
effects of including the NLO fragmentation terms are shown, as well as the effect of the use 
of different charm PDFs on the cross section. In Section IV we summarize our findings. 



II. THEORY 
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FIG. 1: Compton Scattering 



Denoting the electromagnetic and strong couplings by a and a s , respectively, the lowest 
order subprocess for the production of a photon plus a heavy quark is the Compton subpro- 
cess, g + Q —>■ 7 + Q, shown in FigffJ This subprocess is of order aa s and in the variable 
flavor scheme employed here there is only this one subprocess to this order. When one con- 
siders higher order subprocesses such as qQ — > qQ'j, for example, there will be a region of 
phase space where the photon may be emitted collinear with the final state q, giving rise to 
a collinear singularity. This singular contribution may be absorbed into a photon fragmen- 
tation function D 1 / q . The photon fragmentation functions satisfy a set of inhomogeneous 
DGLAP equations, the solutions of which are of order a/a s . More specifically, one has 



dD l/q (z,t) (\>^± 



[D l/q ® P qg + D llg ® P g _ 



dDr f / g (z,t) a 

dt = 2~k [J ~'' '' ' ' iq 7 J/ "' ! '' ' ,/,; 



where t = ln(Q 2 / Aq CD ) and ® denotes a convolution. Writing a s (t) = l/bt it is easy to 
see that the solutions for both D 1 / q and D 7 / g are proportional to both t and a so that the 
fragmentation functions may be thought of as being 0(a/a s ). Therefore, another class of 
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NLO su 


bprocesses 


99 ~ 




gQ- 


igQ 


Qq- 


-> jqQ 


Qq- 


-> jqQ 


qq- 


+ -yQQ 


QQ- 




QQ- 


-* iQQ 



TABLE I: a list of all 2 — > 3 NLO hard-scattering subprocesses 



contributions of order aa s consists of 2 — > 2 QCD subprocesses with at least one heavy 
quark in the final state convoluted with the appropriate photon FF. An example is shown 
in FigfJJ 
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FIG. 2: An example of Leading Order Fragmentation Contributions 1) gg — > where the 

photon can fragment off from either one of the final state heavy quarks, 2) gQ — > where again 
the photon can fragment off from either one of the final state partons, the gluon or the heavy quark 

At the next order in perturbation theory, oioi 2 s , the phase space for producing a photon 
in association with a heavy quark increases and now there are seven possible subprocesses, 
which are listed in Table 1. As in the LO case, there are fragmentation contributions that 
need to be taken into account in order to have a complete NLO calculation. Thus all 2 — >• 3 
QCD subprocesses of order a z s once convoluted with D^/g^^z, Q), give something of the 
NLO: 0(a 3 s ) <g> D 1 / q>g ~ a 3 s a/a s = aa 2 s , FigJSJ 

As soon as we go beyond LO, ultraviolet (UV), soft, and collinear divergences appear. 
The UV singularities arise from virtual diagrams when the momenta of the virtual particles 
go to infinity. To take care of these the theory is renormalized, and the singularities are 
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FIG. 3: An example of Next-to-Leading Order Fragmentation Contributions 1) gg — > QQgj, where 
again the photon is produced by fragmenting from either of the final state partons produced in the 
hard scattering, 2) another example of NLO fragmentation gQ — > ggQ^f 

absorbed into the now renormalized strong coupling a s . The soft singularities appear in the 
case where the energy of a massless particle like the gluon goes to zero, and the collinear 
singularities arise when two massless particles are emitted collinearly. Since generally the 
energies that are considered are much larger than those of the quark masses, the quarks 
are treated as massless, and the calculation is done in the massless approximation. To take 
care of the divergences the calculation is regularized. The regularization scheme that is used 
here is Dimensional Regularization (DR). In DR the scattering amplitudes are computed in 
d = 4 — e dimensions, and the singularities are exposed as poles in 1/e. These poles cancel, 
once the virtual, soft and collinear contributions are added or they have to be absorbed 
into Parton Distribution Functions (PDFs) and FFs with the use of the DGLAP evolution 
equations. Once this is done it is safe to go back to 4 dimensions. 



To perform the NLO calculation the two cutoff phase space slicing method 11] is used. 
In it the phase space is divided into 2^2 and 2^3 body contributions. The 2^3 
body phase space is further divided into a hard region, where no singularities are present, 
a collinear region, where the collinear singularities are present and a soft region, where 
the soft singularities occur. The separation between the different regions is done with the 
help of two parameters : the soft cutoff - 8 S , and the collinear cutoff - 8 C . In the phase 
space slicing method a gluon is considered to be soft, if E g < 5 s y§/2, where \/I> is the 
hard scattering center of mass energy. In order to simplify the integration in this region 
the double pole or soft gluon approximation is used and the 4-momentum of the gluon is 
set to zero, when it appears in the numerator. The collinear region is taken to be where 
either or < S c s, where s^,^- are the Mandelstam variables. In the collinear region 
the leading pole approximation is used, so the relative transverse momentum of the two 
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collinear particles in the numerators of the expansion is neglected. The integration over 
phase space is done with the use of Monte Carlo. This is very useful, as the cross section can 
be calculated as differential in any variable that is needed, such as the transverse momentum 
- px or the rapidity - y, etc. of a given particle, without having to worry about calculating 
different Jacobians. There will be a dependence on the cutoff parameters in both the 2-^2 
contributions (which include the collinear and soft regions) and the 2 — * 3 contributions, 
but this dependence will disappear once the two contributions are added together. 

One final point needs to be addressed concerning the subprocess qq — > jQQ. There is a 
collinear singularity associated with the region where the final state Q and Q are collinear. 
Physically, this corresponds to a ^g final state with the gluon splitting into the QQ pair. 
Normally, this singular region would be integrated over yielding a two-body contribution 
dependent on 5 C which would be proportional to the subprocess qq —>■ jg. This would be 
added to the one-loop corrections for the qq — > jg subprocess, the poles in e would cancel 
and there would be a residual 5 C contribution to the qq — > jg subprocess. This would cancel 
against a similar contribution from qq — ► ^QQ once a suitable jet definition was implemented 
in the calculation. However, once one tags the jet as containing a heavy quark, the problem 
arises in that there is no contribution from the subprocess qq — > jg. Hence, there is an 
uncanceled S c dependence in the 2 — > 3 contribution. This problem is addressed by realizing 
that physically the final state gluon can not fragment into a QQ pair unless its invariant 
mass exceeds 4mQ. Imposing this constraint on the events generated for qq — ► ^QQ avoids 
the problem of the uncanceled 5 C dependence. 

III. RESULTS 

A. Tevatron Predictions 

For the numerical results shown below the CTEQ6.6M PDFs [12j were used, unless oth- 
erwise stated, with a 2-loop a s corresponding to a s (M z ) = 0.118. The cross section was 
calculated for a center of mass energy of v^S* = 1.96 TeV corresponding to the measurements 
being made at the Tevatron. The cuts applied reproduce the ones used by the D0 experi- 
ment, where the lower bounds for the transverse momenta of the photon and heavy quark 
are as follows: py 7 > 30 GeV, ptq > 15 GeV, and their rapidities are limited to the central 
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FIG. 4: The differential cross section, da/dpr-y FIG. 5: contributions of the different subpro- 
for the production of a direct photon and a cesses to the differential cross section, NLO - 

solid line, annihilation qq — > QQ^f - dashed line, 
qQ — ► qQ'y - dotted line, gQ — > gQ'y - dot dashed 
line, gg — > QQ7+LO - dash dot dotted line, 
QQ — * 1QQ1 an< i — * 7QQ - dot dash dotted 
line 

region of the detector \y^\ < 1, < 0.8. If two final state partons lie within a cone of 
radius Ai? = 0.5, where R = \J Arj 2 + A</> 2 then they are merged into a single jet. If a final 
state has two heavy quark jets within the detectable region, it is counted only once, taking 
into account the transverse momentum of the more energetic jet. To be experimentally 
detectable, a photon needs to be isolated. This means that it should not be surrounded by 
hadronic energy more then — e * E 1 in a cone of radius R = Ri SO around it. The photon 
isolation requirements imposed model those needed in the DO detection of a direct photon 
and are: R x < 0.2, e 1 < 0.04 and R 2 < 0.4, e 2 < 0.07. 

The differential cross section for the process pp — > 76X as a function of the transverse 
momentum of the photon is shown in FigJH It is interesting to note in FigJH that as pr 7 
grows, the difference between the LO and NLO curves increases substantially. To understand 
the origin of this effect it is necessary to look at how the different subprocesses listed in Table 
1 contribute to the cross section. This decomposition is shown in FigHJ 

It is apparent from Figj5] that the effect shown in FigJH is driven by the annihilation 
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subprocess, qq — > 'yQQ, which overtakes the Compton contribution and starts dominating 
the cross section at pr-y ~ 70 GeV. The gQ — > jgQ and Qq — > jqQ /Qq — > jqQ subprocesses 
contribute to the cross section about equally, with gQ — > jgQ prevailing over Qq — > jqQ 
I Qq 1Q.Q a t small pt^, where the gluon PDF is larger than the light quark PDF, and then 
at large Pt^, Qq ^qQ /Qq iqQ takes over when the light quark PDFs become larger 
than the gluon PDFs. The Compton subprocess and gg — > •yQQ are added together, since 
the gg — ► ^QQ contribution is negative. This negative contribution is what remains after 
the appropriate collinear terms are subtracted. The size of the 2 — > 3 gg contribution is scale 
dependent, with the compensating collinear terms contributing to the 2 — > 2 component. 
The role of the QQ — > •yQQ / QQ — > jQQ subprocesses is almost negligible, as the heavy 
quark PDFs are much smaller than the light quark and gluon PDFs. 

Since the annihilation subprocess dominates the cross section at large pr-y, it is useful 
to look at some of the Feynman diagrams contributing to it, as shown in Figj6j There 




FIG. 6: Some typical Feynman diagrams for the annihilation subprocess qq — ► "yQQ where 1) the 
photon is emitted from the final state heavy quarks and 2) the photon is emitted from the initial 
state light quarks 

are two channels through which the annihilation subprocess can be produced, an s-channel 
shown in diagram 1) and a t-channel in diagram 2) of FiglBJ Since the photon couples to 
the final state heavy quarks in diagram 1), this diagram is proportional to the heavy quark 
charge, eg, whereas in diagram 2) the photon couples to the initial state light quarks, and 
thus this diagram is proportional to the light quark's charge, e q . Diagram 2) begins to 
dominate as px^ grows and, since it does not depend on the heavy quark charge, we expect 
the difference between the bottom and charm cross sections to decrease as pt 1 increases. 
This indeed is the case as can be seen from Figf?l where the NLO differential cross section 
for the charm quark (solid line) and the one for the bottom quark (dashed line) tend to come 
closer as the value of the transverse momentum increases. However, the difference between 
the LO cross sections stays about the same as can be seen from the dot-dashed (charm) and 
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FIG. 7: a comparison between the differential FIG. 8: the ratio of the charm and bottom differ- 

cross sections, da/dp?^ for the production of a ential cross sections versus pt^, at NLO - solid 

direct photon and a bottom quark, and that of line and at LO - dashed line 

a direct photon plus a charm quark at NLO and 

LO, charm at NLO - solid line and for bottom 

at NLO- dashed line, charm at LO - dot dashed 

line, bottom at LO - dotted line 

dotted curves (bottom) in FigJTJ, and also from FiglHl where the ratio of the NLO and LO 
charm and bottom differential cross sections is shown. The ratio of the LO cross sections 
stays almost constant since the main contribution to the LO cross section comes from the 
Compton subprocess, with the difference between the charm and bottom curves arising from 
the difference in the charges of the charm and bottom quarks and the relative sizes of the 
heavy quark PDFs. The ratio of the two LO cross sections depends on the ratio of the 
charges squared which is e 2 c /el = 4, and is driven up from that value to about ~ 7 because 
the charm PDF is larger than the bottom PDF. 

The fact that the annihilation subprocess dominates the cross section at large pt^ also 
increases the scale dependence of the cross-section in that region. There is no Born term 
which involves a qq initial state and, therefore, the contributions from the annihilation 
subprocess start in 0(aa 2 s ). As such, the typical compensation between LO and NLO 
contributions for this subprocess is missing and the annihilation subprocess can be thought 
of a leading order. 
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FIG. 9: Scale dependence of the differential cross section, da/dpT-y for the production of a direct 
photon and a bottom quark, where the three different scales have been set to be equal \i = \i r = 
[if = fj,p, fi = px-y - solid line, /i = pt~//2 - dashed line, fj, = 2px-y - dotted line 

As can be seen from Fig J9] the scale dependence increases at large pr 7 , where the annihi- 
lation starts to dominate. The renormalization, fi r , factorization, fif and fragmentation, fip 
scales have been set to be equal and are denoted by /i. 

B. LHC Predictions 

When data become available from the Large Hadron Collider (LHC), it will be very 
important to have a good understanding of what Standard Model (SM) processes are going 
to look like at center of mass energies of = 14 TeV, as these processes will provide 
important means of calibrating and understanding the detectors and, ultimately, are likely 
to provide significant backgrounds to new physics signals. The differential cross section 
versus the transverse momentum of the photon is shown in FigJTOl It is apparent that the 
increase of the difference between the NLO (solid line) and the LO (dashed line) grows much 
less rapidly with increasing px-y than was the case for the Tevatron. In FigfrT] the K-factor, 
which is the ratio between the NLO and LO cross section for b quarks is shown, which 
stays stable and is around 2. To understand the difference between the LHC and Tevatron 
curves, the contributions of the different parts contributing to the LHC cross section are 
shown in Figfl2l From Figfl2] it can be seen that the annihilation subprocess no longer 
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FIG. 10: the differential cross section versus the FIG. 11: K factor, or the ratio of the NLO to 

transverse momentum of the photon da/dpT^ the LO differential cross section for pp — > bjX 

for the production of a direct photon and a at = 14 TeV 

bottom quark at LHC center of mass energies, 

y/S = 14 TeV, NLO - solid line, LO - dashed 

line 

drives the cross section at high pxy, and now it is the LO, and the gQ — > 'jgQ subprocesses 
that are the most prominent. These differences come about for two reasons. As the LHC 
collides two beams of protons, instead of the proton and antiproton beams at Fermilab, 
there is no longer any valence light antiquarks present. Hence, the relative contribution of 
the annihilation subprocess is decreased. Also, because the LHC will ultimately operate at 
a center of mass energy which is about seven times larger than that of the Tevatron, lower 
values of x ~ Pr/y/s are probed at the LHC. For the kinematic region shown in Fig. [12] 
the gluon PDF is dominant, accounting for the continued importance of the gQ initiated 
subprocesses. 

An interesting consequence of this pattern of subprocess contributions is that the the 
dominant parts are all proportional to the heavy quark PDFs. Such was not the case for the 
Tevatron curves, except for the low end of the pt-j range. Accordingly, heavy quark + photon 
measurements at the LHC will have the potential to provide important cross checks on the 
perturbatively calculated heavy quark PDFs. These PDFs are likely to provide important 
contributions to other physics signals - either standard model or new physics - and such 
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FIG. 12: contributions of the different subprocesses to the differential cross section, da /dpx^ for 
pp — ► bjX at y/S = 14 TeV, NLO -solid line, annihilation qq — > QQj - dashed line, qQ — > qQ^y 
- dotted line, gQ — > gQ^ - dot dashed line, gg — ► QQ7+LO - dash dot dotted line, ~~ * iQQ, 
and — ► 'yQQ - dot dash dotted line 

checks will be an important part of the search for new physics. 
C. NLO Fragmentation and Photon Isolation 

It is interesting to investigate what effect the NLO Fragmentation contributions have 
upon the cross section. FigIT3l shows the ratio between the full NLO calculation and the 
cross section with only LO fragmentation. If there are no isolation requirements imposed 
on the photon, the cross section increases up to ~ 30%, solid curve in FigfTBl As mentioned 
above a photon needs to be isolated in order to give a clear signal at a detector. The isolation 
requirements affect the photon which is produced by fragmentation the strongest, as it is 
emitted in close proximity to the parton from which it is fragmenting. This can be seen from 
the dashed line in Fig. [13], where the NLO fragmentation contribution has now decreased 
to a few %. FigJTJ] shows the comparison between the differential cross section with the 
inclusion of isolation and without it. As can be seen this difference is larger at low pr 7 , 
but the two curves come close to each other with the increase of the photon's transverse 
momentum, where as seen from Fig J5] the qq — > 7QQ subprocess takes over the cross section. 
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FIG. 13: ratio between the differential cross sec- 
tion da/dpT'y, with NLO fragmentation contri- 
bution included and the differential cross section 
with just LO fragmentation included, solid line 
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FIG. 14: comparison between the differential 
cross section, da/dpT-y without isolation require- 
ments and with them, no isolation - solid line , 
isolation - dashed line 



In the CTEQ6.6M PDFs used in the previous sections, the charm quark is radiatively 
generated from the gluon's PDF with the use of the DGLAP equations. Thus it follows that 
there are no charm quarks present at scales below the charm mass,m c , or that the charm 
PDF, c(x,fi) = 0, when /i < m c . This however does not need to be the case, and there are 
models that study the possibility for an intrinsic charm component of the nucleon [3]. Two 
such models are the BHPS model, which is a light-cone model, and the sea-like model in 
which the charm distribution follows the shape of the light flavor sea quarks. The difference 
between the three cases is shown in FigJT5| where the solid curve shows the CTEQ6.6M or 
radiatively generated charm scenario, the dashed curve is the CTEQ6.6C2, or BHPS model, 
and the dotted curve is CTEQ6.6C4 PDF or the sea-like model. The difference between the 
BHPS distribution and the radiatively generated case are most noticeable at large x, whereas 
the sea-like model is about equally larger than the CTEQ6.6M PDF at all values of x. How 
these different PDFs affect the cross section can be seen from FigJT6l The dotted curve 
shows the cross section generated with the use of the sea-like intrinsic charm PDF, and it 
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FIG. 16: the differential cross section, da/dpx-y, 
for the production of a direct photon and a 
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is larger than the solid curve by about the same amount at all values of pr 7 . The difference 
between the radiatively generated charm cross section and the BHPS charm however is not 
great at small transverse momentum, but it increases at large PT-y, as is expected given the 
differences between the CTEQ6.6M and CTEQ6.6C2 PDFs at large x. 



IV. CONCLUSION 



We have presented the results for the inclusive cross section for the production of a 
direct photon in association with a heavy quark, pp/pp — > jQX, up to 0(aa 2 s ) with NLO 
fragmentation included. The inclusion of NLO fragmentation has a noticeable effect on 
the cross section if no isolation is imposed. However, this effect decreases if isolation cuts, 
needed for a clean photon signal, are imposed. Predictions were presented for pp collisions 
at v^S* = 1.9QTeV and for pp collisions at 14 TeV. At the Tevatron, due to the pp beams, 
the valence quarks and antiquarks are dominant, and thus it is the annihilation subprocess 
qq — ► jQQ that dominates the cross section at large pt t Therefore the sensitivity to the 
initial state heavy quarks and their content in the nucleon decreases, and the difference 
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between the bottom and charm differential cross sections, da/dpx-y, also diminishes. At the 
LHC, where two beams of protons are colliding, there are no longer any valence antiquarks 
present, and processes with initial gluons and heavy quarks dominate. Thus there should be 
a greater possibility to learn more about the heavy quark role in the nucleon. In particular, 
the perturbatively calculated heavy quark PDFs may be checked using such data. 
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